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Current algorithms used to put a lattice gauge configuration into Landau gauge either suffer from the problem 
of critical slowing-down or involve an additional computational expense to overcome it. Evolutionary Algorithms 
(EAs), which have been widely applied to other global optimisation problems, may be of use in gauge fixing. Also, 
being global, they should not suffer from critical slowing-down as do local gradient based algorithms. We apply 
EA's and also a Steepest Descent (SD) based method to the problem of Landau Gauge Fixing and compare their 
performance. 
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1. Motivation 

As SD is a local algorithm the time taken to 
converge goes as some power (> 1) of V the lat- 
tice volume (critical slowing down) 0] . A solution 
to the problem is Fourier acceleration but the 
cost of this goes as the cost of a FFT - VlogV at 
best. Once the gauge fixing condition has been 
met, Gribov copies may still exist. These can con- 
tribute to numerical errors in the measurement of 
various gauge dependent quantities, || @ [|) @ 
most notably propagators. 

EA's on the other hand are global algorithms 
and should not suffer from critical slowing down. 
They are capable of finding a global maximum 
and so by suitable choice of constraint they can 
in principle eliminate Gribov copies. 

2. Landau-gauge fixing on a lattice 

For Landau-gauge we seek 



= 



(1) 



A sufficient (but not necessary) condition for this 
to be true is the maximisation of 



F 



afxA a (x)A sll {x) 



where 

Al(x) = A lt (x)-d* X (x) 



(2) 
(3) 



Following |2| we do the following: 

• Name the lattice version of d^A^, O, and 
the lattice version of F, F^. 

• Implement SD. 

• Compare the performance of SD with a 
competing algorithm in maximising Fl , us- 
ing 6 as a measure of convergence. 

3. EA's applied to Landau-gauge fixing on 
a lattice 

EAs solve global optimisation problems with 
natural selection and genetic style operations. 
The idea is to evolve successive generations of 
genomes with the goal of maximising their fitness 
01 @ ||- We now describe the implementation 
of these operations to gauge fixing. 

3.1. Coding 

We encode the gauge transformations G(x) 
(not the gauge fields U(x)). Each genome is made 
of group elements for each site of the lattice (ID is 
shown in Figure [j]) . Each group element is made 
of floats which represent the angles parameteris- 
ing it. 

• U(l) - just one angle. 

• SU(2) - use the angles in 0(3). 
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Figure 1. Encoding of groups 
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Figure 2. Evaluating the fitness of a gauge trans- 
formation 



• SU(3) factorise into 3 SU(2) matrices acting 
on the three SU(2) subspaces in SU(3). 

3.2. Mutation 

The floats are mutated by adding some nor- 
mally distributed random numbers to them. The 
standard deviation of the distribution can vary 
according to the generation, genome and position 
within the genome. 

3.3. Recombination 

Cutting and splicing of the genome as used in 
Ref. can be used. Best results are achieved by 
cutting on group boundaries and by having the 
genome duplicate the dimensionality and bound- 
ary conditions of the lattice. More successful 
methods are those which globally interpolate be- 
tween genomes. 

3.4. Fitness 

The fitness of a genome is found by applying 
the G(x) that it represents to the original gauge 
fields and then evaluating Fl on the result as 
shown in Figure |[ 

Genomes are selected from a population for 
breeding using algorithms which favour those 
with high fitness while at the same time main- 
taining diversity in the population. The last con- 
dition is to allow the population to escape local 
maxima. 

4. Typical EA results 

• The graphs that follow are not measures of 
flops. 



• One generation of EA takes more time than 
one iteration of SD. 

• There is also the issue of memory. For an 
EA one must store many genomes, each of 
which has the size of a gauge transformation 
for the entire lattice. 

• All lattices used are 8 1 with U(l) gauge 
fields with the exception of the last two fig- 
ures which are 16 1 with U(l) gauge fields. 

Figure |^ shows the statistics of a population as 
the algorithm progresses. The standard deviation 
of the fitness gives a measure of convergence (in 
addition to 0). Figure ^| shows a direct compari- 
son between SD and EA. On average SD performs 
better, and we have always found this to be the 
case. 

The reason for this might be that there is no 
gradient information available to EA's. Close 
to convergence the genomes are doing a random 
walk (ie. without guidance) in the phase space 
of gauge transformations. Surprisingly, spectral 
decomposition of the gauged configurations sug- 
gests that the EA's are suffering from a similar 
problem that afflicts SD {§ . Examination of Fig- 
ures H and H shows that the slow decay of low 
momenta modes in both algorithms is a problem. 

5. Conclusion 

In this preliminary study, our efforts have con- 
centrated on improving the rate of convergence 
of EAs. This seems to go as some power of r/n 
where r is the distance from the actual maxi- 
mum and n is the number of floats in a genome. 
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Figure 3. Fitness of Population vs Generations 
a typical plot for EA 
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Figure 4. 6 vs Generations - a comparison of 
the two methods for a typical run 
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Figure 5. Power vs Iterations - Fourier compo- 
nents of successive gauged fields using SD 



Figure 6. Power vs Generations - Fourier com- 
ponents of successive gauged fields using EA 



Presently we are looking at their application to 
distinguishing between Gribov copies as in the- 
ory EAs should be able to find the global max- 
imum effectively. Various hybrid methods com- 
bining EAs and local gradient methods - such as 
SD or its variations - are also being considered. 
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